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Estimates are made of ultra-high energy neutrino cross sections based on an extrapolation to 
very small Bjorken x of the logarithmic Froissart dependence in x shown previously to provide an 
excellent fit to the measured proton structure function F%(x, Q 2 ) over a broad range of the virtuality 
Q 2 . Expressions are obtained for both the neutral current and the charged current cross sections. 
Comparison with an extrapolation based on perturbative QCD shows good agreement for energies 
where both fit data, but our rates are as much as a factor of 10 smaller for neutrino energies above 
10 9 GeV, with important implications for experiments searching for extra-galactic neutrinos. 
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Introduction: The experimental effort to detect extra- 
galactic, ultra-high-energy (UHE) neutrinos has grown 
rapidly in the past decade. Optical [l[ and radio 
telescopes and cosmic ray air shower arrays 0] are now 
searching for evidence of point and diffuse neutrino 
sources up to and beyond EeV energies. Proposals have 
been made and others are in preparation Q for new tele- 
scopes or expansions of ones currently deployed, and am- 
bitious satellite-born telescopes have been proposed Q. 
The highest energies proposed reach beyond 10 12 GeV. 

Critical to all of this effort are accurate estimates of 
event rates, based on the extrapolation of measured neu- 
trino deep-inelastic scattering (DIS) cross sections to en- 
ergies far beyond currently available data 0,0,1- The 
estimates are only as reliable as the extrapolations, and 
determination of fluxes and extraction of signals of new 
physics at UHE depend on them. Most existing extrap- 
olations are done within the framework of perturbative 
quantum chromodynamics (pQCD), and they involve ex- 
tending fitted parton distribution functions (PDFs) into 
domains in Bjorken x much below those now accessible 
experimentally, and into domains in which linear pQCD 
evolution [9j is of questionable applicability. Other phys- 
ical phenomena are expected to alter the x dependence 
in this very small x region [loj ]. although a complete an- 
alytic solution does not yet exist. 

New, alternative methods of extrapolation in x are 
of significant interest, both theoretically and for phe- 
nomenological applications. Imposition of the Frois- 
sart II l( unitarity and analyticity constraints on inclusive 



deep-inelastic cross sections [12] leads to the expectation 
that the x dependence of the proton structure function 
F^XjQ 2 ) should grow no more rapidly at very small x 
than ln 2 (l/a;). This relatively slow growth may be con- 
trasted with the more rapid inverse power dependence 
characteristic of PDFs. Excellent fits to data were ob- 
tained for x < 0.1 with an assumed logarithmic ex- 
pansion, for a wide range of virtuality Q 2 . We explore 



in this Letter the consequences of the Froissart logarith- 
mic form for UHE neutrino phenomena, computing both 
neutral and charged current cross sections. In doing so, 
rather than working with parton distribution functions 
for the decomposition into quark and antiquark contri- 
butions, we devise and test a procedure based directly on 
experimental F% data. We obtain excellent agreement 
with extrapolations based on the CTEQ4-DIS parton 
densities in the neutrino energy range less than 10 8 GeV. 
However, we predict an important departure for larger 
energies, with our neutrino cross sections being about a 
decade smaller at the highest energies. At the very least, 
our results suggest that estimates that fall between ours 
and those obtained from PDF extrapolations be used for 
guidance in the consideration of new experiments. 

Neutrino-isoscalar nucleon cross sections: In the stan- 
dard parton model the inclusive differential cross section 
for the charged current (CC) reaction vi + N — > l~ + X 
on an isoscalar nucleon N = (n + p)/2 and the neutral 
current (NC) cross section i/g + N —* vi + X , where in 
both cases, I = e, ji, t, is 
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where — Q 2 is the invariant squared momentum transfer 
between the incoming neutrino and the outgoing muon, 
m is the proton mass, and Gf is the Fermi coupling 
constant. The intermediate vector boson mass, My, is 
M w = 80.4 GeV for CC and M z = 91.2 GeV for NC. 
Symbols qi and i =CC, NC, are linear combinations 
of quark and antiquark PDFs. The Bjorken scaling vari- 
ables, where v = E v — Eg is the energy loss in the labo- 
ratory frame, are given by 



2mv ' 



y 



V 

E~u 



< x,y < I. 



(2) 



Charged current cross section: With valence and sea 
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quark distributions denoted by subscripts v and s, re- 
spectively, the relevant PDFs in Eq. JT]) are 



qcc(x,Q ) 



u v (x,Q 2 ) +d v (x,Q 2 ) 



u s (x,Q 2 ) + d s (x,Q 2 ) 
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and 



qcc{x, Q ) 
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+c s (x,Q 2 )+t s (x,Q 2 ), (4) 



where u, d, c, s, i, and b represent the contributions from 
the up, down, charm, strange, top, and bottom flavors. 

Neutral current cross section: The relevant PDFs in 
Eq. |T]) involve chiral couplings L u — 1 — | sin 2 6w , Ld = 
— 1 + | sin 2 9w, R u = — § sin 2 i?^ = | sin 2 6\y, where 
sin 2 8w = 0.226 is the weak mixing parameter. For de- 
tails, see Ref. [7]. 

Kinematics: Replacing Q 2 in Eq. ([1]) by Q 2 = 
2mE v xy, we obtain an expression in terms of E v ,x and 
y. We choose to integrate first over y. To avoid singular- 
ities in the integration, we introduce Q 2 lin = 0.01 GeV 2 , 
such that Q 2 = 2mE u xy > Q 2 nin - This defines x min , 
the x-integration minimum, as x min = Q 2 ain /(2mE v ). 
Thus, for Xmin < x < 1, our integration limits for y 
are y min = x min /x < y < 1. 

The vector boson propagator, (M V /{Q 2 + M v )) , es- 
sentially fixes an "effective" x at x c h ~ M v /(2mE„). For 
E v = 10 12 GeV, this means we must explore quark dis- 
tributions having x c s ~5x 10 -9 , at Q 2 ~ M v ~ 10, 000 
GeV 2 , both of which involve enormous extrapolations 
from currently available structure function data. At these 
energies, the propagator also serves to make the calcula- 
tion insensitive to the choice of Q^^. 

Analytic expression for the structure function: In prior 
work [12j, it was shown that an excellent fit to the DIS 
structure function for x < xp, is given by 
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where 



A(Q 2 ) = a + a x In Q 2 + a 2 In 2 Q 2 
B(Q 2 ) = bo + hlnQ 2 + 6 2 ln 2 Q 2 . 



(6) 



The fitted numerical values of a.j and bk and their un- 



certainties may be found in Ref. [12j; Fp = 0.41, and 
2- P = 0.09. 

The bulk of the neutrino cross section comes from ex- 
ceedingly small x. For large x, where xp < x < 1, it 



suffices to approximate the proton structure function by 



Fi(x,Q 2 ) = 



F P 



^ 2 \l~xf, (7) 



where the exponent a(Q 2 ) is chosen so that the first 
derivatives of Eq. © and Eq. (O are equal at x — xp. 
This choice satisfies the spectator valence quark count- 
ing rule [l3| F^ix) — > as (1 — a;) 3 as x — > 1. Numerical 
analysis shows that this choice has the important conse- 
quence that the integral of the proton structure function 
over x is nearly constant over an enormous Q 2 range, i.e., 



f F%(x,Q 2 )dx w0.16, 0.1 < Q 2 <10 5 GeV 2 
Jo 



(8) 



The constant 0.16 is compatible with results that show 
that quarks carry ~ 50% of the momentum in a proton. 

The description of F 2 p (x,Q 2 ) by Eqs. (5) - (8) yields 
a high quality fit to the HERA inclusive deep-inelastic 
data for all x and Q 2 . 

"Wee parton" picture: We obtain the quark distribu- 
tion functions in Eq. ([T]) from a wee parton model for very 
small Bjorken x, having the following features: 

• there are essentially only sea quarks at small 
enough x, with negligible valence quark contribu- 
tions (for earlier use, see Ref. Q), i.e., we set 

u v (x,Q 2 ) = d v (x,Q 2 ) = 0. 

• all sea quarks give equal contribution (i.e., equipar- 
tition), 

U(x,Q 2 ) = u s {x,Q 2 ) = u s (x,Q 2 ) = d s (x,Q 2 ) = 
d s (x,Q 2 ) = s s (x,Q 2 ) = s s (x,Q 2 ) = c s (x,Q 2 ) = 
c s (x, Q 2 ). 

If only two families contribute (u, d, c, and s), 
F£(x, Q 2 ) =V] e 2 x[ gi (x, Q 2 )+g 4 (x, Q 2 )],i = 1, ... 4, (9) 



or, alternatively, 



xU(x,Q 2 ) = —F£(x,Q 2 



(10) 



for x < x max , where x max ~ 10~ 3 — 10~ 4 . If we had 
used only one family of quarks — u, d — or three families — 
u, d, c, s, t, b — instead two families — u, d, c, s — we would 
also find that xq(x,Q 2 ) = xq(x,Q 2 ) = ^F%(x, Q 2 ), so 
that Eq. (Q]) for charged currents is independent of the 
number of families. A similar result is true for the neu- 
tral current cross section. Employing this picture, we 
find that accurate knowledge of F^ix^Q 2 ) at small x 
and large Q 2 provides the ingredients necessary to calcu- 
late the charged and neutral current neutrino cross sec- 
tions. The fitted form of Eq. §5§ is sufficiently accurate to 
furnish us with quark distribution functions having the 
needed precision. Using the full squared error matrix for 
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FIG. 1: Charged (CC) and neutral (NC) current neutrino 
cross sections in cm 2 vs E„, the neutrino energy in GeV. The 
solid and dash-dot-dot curves are our CC and NC cross sec- 
tions, respectively, for 10 < E v < 10 14 GeV, based on a 
proton structure function that varies as ln 2 (l/:r) for small x. 
The long dash curve and the dash-dash-dash curve are the 
Gandhi et al. 0] CC and NC cross sections , respectively, 
for 10 < E v <10 12 GeV, based on the CTEQ4-DIS quark 
distributions. 



= 4(1 - 2 sin 2 (9™ + — sin 4 6 w )xU(x,Q 2 ) 
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= 2MxU(x,Q 2 ) = l.l9F£(x,Q 2 ), (12) 

where Eq. (fTtJ]) is used in the last step. The neutral cur- 
rent cross section simplifies considerably. For direct com- 
parison with the charged current cross section of Eq. (Ill)) , 
it can be rewritten as 

c^onc v _ 2& F mE v ( Mf V 
dxdy [ v) ix \Q 2 + M 2 ) X 

[0.298if(z, Q 2 )} (2-2y + y 2 ). (13) 

To the extent that the Z propagator is somewhat less 
restrictive as a cutoff than the W propagator, comparison 
of Eq. (T3J) and Eq. ((TTJ) shows that the ratio of the NC 
cross section to the CC cross section is > 0.298/0.9 = 
0.33, independent of energy. Numerical evaluation gives 
0.40 at E v = 10 7 GeV, slightly higher because of the Z 
propagator. Our NC cross section for isoscalar nucleons 
is given in Table Q] and shown in Fig. [T]as the dash-dot- 
dot curve, plotted in the energy interval 10 < E v < 10 14 
GeV. The Gandhi et al. [7] NC cross section, for 10 < 
E v < 10 12 GeV, is the dash-dash-dash curve. Again, the 
agreement is excellent up to E v ~ 10 s GeV. 



the structure function determination [12|, we find that 
F%(x = 10" 8 ,g 2 = 6400 GeV 2 ) = 24.84 ± 0.17, a frac- 
tional statistical accuracy of only ~ 0.7%. This very small 
uncertainty due to parameter errors assumes, of course, 
the validity of our ln 2 (l/a;) model at very small x. 

Charged current cross section evaluation: For our 
model, xq C c{x,Q 2 ) = xq cc {x,Q 2 ) = 2xlI(x,Q 2 ). Thus 
xU(x,Q 2 ) = (z, Q 2 ) and Eq. P simplifies to 
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with F^x^Q 2 ) given by Eq. ([5]) for < x < xp and 
Eq. for x P < x < 1. 

Results of a direct double integration of Eq. (fTTj) . with 
Qmin = 0-01 GeV 2 , for the neutrino energy range 10 < 
E v < 10 14 GeV, are given in Table |T] and shown in Fig. 
[T] as the solid curve. Also shown, for comparison, are 
the results of Gandhi et al. [7] for the CC cross section 



with the quark distributions from CTEQ4-DIS 14]. The 
Gandhi et al. curve — the long dash curve — covers the 
energy range 10 < E v < 10 12 GeV. The agreement up to 
neutrino energies <, 10 s GeV is striking. 

Neutral current cross section evaluation: For our 
model, the NC quark distributions in Eq. ([T]) are 



xq NC (x,Q 



xq NC (x,Q 2 ) = 



2xU(x,Q 2 ) x 



TABLE I: Neutrino CC and NC total cross sections, with 
neutrino energy E v energy in GeV and cross sections in cm 2 . 
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Robustness of cross sections: The differential cross sec- 
tions were evaluated numerically in Mathcmatica and 
found to be numerically stable, essentially independent of 
Q^ijj and the methods of integration. The dependence of 
the cross sections on the functional form of F^iQ 2 ,x) for 
1 > x > xp was tested by setting Ff {Q 2 , x) ~ x(l — x) 3 
for large x, and the change was found to be ~ 2% 
at E u = 10 s and ~ at E v = 10 12 GeV. If we set 
F2(Q 2 ,x) = for 1 > x > xp, an extreme case, we 
find the changes to be 6% at E u — 10 8 GeV and ~ 
at E v = 10 12 GeV. We tested our equipartition hypoth- 
esis by changing the strengths of the heavy sea quark 
distributions such that 

s s (x,Q 2 ) = s s {x,Q 2 ) = 0.96 U(x,Q 2 ) 
c s (x,Q 2 ) = c s (x,Q 2 ) = 0.80 U(x,Q 2 ), (14) 

similar to the distributions used by CTEQ. This change 
gives us cross sections that are ~ 6% greater at E v = 10 8 
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GeV and - 3% greater at E v = 10 12 GeV. These vari- 
ations are negligible compared to the very large differ- 
ences with respect to the cross sections of Gandhi et al. 
Q at the highest neutrino energies. Our calculations 
are numerically stable with regard to our choice of x m - m 
in the integration, and thus, insensitive to our choice of 
Q min = 0.01 GeV 2 

Conclusions: We compute ultra-high energy neutrino 
cross sections based on an extrapolation to very small 
Bjorken x of the logarithmic Froissart dependence in x 
shown previously to provide an excellent fit to the mea- 
sured proton structure function F^x.Q 2 ) over a broad 
range of the virtuality Q 2 . In order to devise expressions 
for the neutral current and the charged current cross sec- 
tions, we first extract quark and antiquark contributions 
based on a simple equipartition wee parton picture valid 
for X max < 1Q- 3 - 10~ 4 or E™ in > 3 x 10 6 - 3 x 10 7 GeV. 
However, it is gratifying to see in Fig. [T]that we are in 
excellent agreement with calculations based on CTEQ4- 
DIS parton densities over the much larger energy range 
10 < E v < 10 s GeV. The two sets of expectations di- 
verge for E v > 10 s GeV, as may be expected since our 
proton structure functions agree with those from CTEQ 
only for x- values greater than 10~ 3 [12j|. The increasing 
differences for x < 10~ 3 reflect the fundamental differ- 
ence in the assumed functional forms for the x depen- 
dence, in our case a form that is constrained to increase 
no more rapidly than ln 2 (l/:r), in contrast to the inverse 
power growth in the CTEQ case. For large neutrino 
energies — above 10 9 GeV — where much smaller x is sam- 
pled, our Froissart-bound-model neutrino cross sections 
are as much as a decade smaller than those based on a 
pQCD extrapolation, a consequence of the fact that our 
structure function F^x^Q 2 ) is significantly smaller at 
small x. The very small x region is also the region where 
our wee parton picture is most robust. 

The region of very small a; is a region of growing inter- 
est theoretically It is a region in which non-perturbative 
physics is expected to set in [l(| and in which linear 
DGLAP pQCD evolution is not expected to hold. While 
we cannot claim that logarithmic dependence on x will 
result from a first-principles solution to small x dynamics, 
neither can we expect an inverse power form to survive. 
The logarithmic form we use provides an excellent fit to 
data over the range in x and Q 2 where it has been tested. 
Its extrapolation to energies relevant in UHE neutrino 
studies provides estimates for event rates that should be 
taken into serious consideration for the planning and data 
analysis of new experiments. 

Acknowledgments: We thank the Aspen Center for 
Physics for hospitality during the research for this 
manuscript. The U. S. Department of Energy, Division 
of High Energy Physics, under Contract No. DE-AC02- 
06CH11357 supports E.L.B. The U. S. Department of 
Energy under Contract DE-FG02-91ER40688, TASK A 



supports in part C-I.T. D.W.M. is supported by DOE 
Grant No. DE-FG02-04ER41308. M.M.B thanks Prof. 
A. Vainstein of the University of Minnesota and Prof. L. 
Durand of the University of Wisconsin for both lengthy 
and fruitful discussions and encouragement. 



[1] A. Achterberg et al. [IceCube Collaboration], Phys. Rev. 
D 76: 027101 (2007), 507 (2004); A. Achterberg et al. 
[IceCube Collaboration], arXiv:astro-ph/0705.1315; G. 
Anassontzis et al. [NESTOR Collaboration], Nucl. Inst. 
Meth.479 439 (2002); J. A. Aguilar et al. [Antares Col- 
laboration], Nucl Inst. Met. A 570, 107 (2007); C. Spier- 
ing et al, [Bailkal collaboration], Nucl. Phys. B Proc. 
Suppl. 138, 175 (2005). 

[2] I. Kravchenko et al. [RICE Collaboration], Phys. Rev. D 
73, 082002 (2006); S.W. Barwick et al. [ANITA Collab- 
oration, Phys. Rev. Lett. 96, 171101 (2006); N. Letinen 
et al. [FORTE Collaboration], Phys. Rev. D 69, 013008 
(2004); P. Gorham et al, [GLUE Collaboration], Phys. 
Rev. Lett. 93, 041101 (2004). 

[3] "Limits on the Diffuse Neutrino Flux from the HiRes 
Data", W. Deng et al [HiRes Collaboration], Proceed- 
ings of the 29th International Cosmic Ray Conference; 
X. Bertou et al, Astropart. Phys. 17, 183 (2002). 

[4] P. Gorham et al., Phys. Rev. D 72, 23002 (2005); H. 
Landsman [AURA Collaboration], Nucl. Phys. B Pro- 
ceedings Suppl. 168, 268 (2007) 

[5] S. Bottai [EUSO Collaboration], Proceedings of the 27th 
International Cosmic Ray Conference, Aug. 2001, Ham- 
burg, p. 848 (2001) 

[6] Yu. Andreev, v. Berezinsky and A. Smirnov, Phys. Lett. 
B 84, 247 (1979); M. H, Reno and C. Quigg, Phys. Rev. 
D 37, 657 (1987); R. Gandhi, C. Quigg, M. H. Reno and 
I. Sarcevic, Astropart. Phys. 5, 81 (1996). 

[7] R. Gandhi, C. Quigg, M. H. Reno and I. Sarcevic, Phys. 
Rev. D 58, 093009 (1998). 

[8] D. W. McKay and J. P. Ralston, Phys. Lett. B, 187, 103, 
(1986); G. M. Frichter, D.W. McKay and J. P. Ralston, 
Phys. Rev. Lett. 74, 1508, (1994). 

[9] V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 
438 (1972), 15, 675 (1972); Yu. L. Dokshitzer, 
Sov. Phys. JETP 46, 641 (1977); G. Altarelli and 
G. Parisi, Nucl. Phys. B 126, 298 (1977). 
[10] L.V. Gribov et al, Phys. Rep. 100, 1 (1983); A. Capella 
et al, Phys. Rev. D 63, 054010 (2001); G. Soyez, 
Ph ys. Rev. D 71 076001 (2005); R. C. Brower et 
al, |hep-th/0603115| larXiv:0707.2408l [he p-th] ; Y. Hatta, 
E. Iancu, and A. H. Mueller, larXiv:0710.2T48l [hep-th] . 
[11] M. Froissart, Phys. Rev. 123, 1053 (1961). 
[12] M. M. Block, E. L. Berger, and C-I Tan, Phys. Rev. 
Lett. 97, 252003 (2006); E. L. Berger, M. M. Block, and 
C-I Tan, Phys. Rev. Lett. 98, 242001 (2007). 
[13] R. Blankenbecler and S.J. Brodsky, Phys. Rev. D 10, 
2973 (1974); J.F. Gunion, Phys. Rev. D 10, 242 (1974); 
S.J. Brodsky and G.P. Lepage, in Proc. 1979 Summer 
Inst, on Particle Physics, SLAC (1979). 
[14] CTEQ Coll aboration, H. Lai et al ., Phys. Rev. D 55, 
1 280 ( 1997) . |arXiv:hep-ph/0611254| 



